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We study a number of hierarchical network models related to the Chalker-Coddington model of 
quantum percolation. Our aim is to describe the physics of the quantum Hall transition. The 
hierarchical network models are constructed by combining series and parallel composition of quantum 
resistors. The localization-delocalization transition occurring in these models is treated by real 
space renormalization techniques. Essentially, the localization-delocalization transition is due to a 
competition between two one-dimensional localization mechanisms. 

PACS numbers: 73.40.Hm, 64.60.Ak 
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I. INTRODUCTION 

Despite the explosion of interest in and seemingly in- 
exhaustible richness of the quantum Hall effect (QHE), 
surprisingly little progress has been made on what is ar- 
guably the central phenomenon associated with all in- 
compressible quantum Hall liquids - the quantum Hall 
transition itself. By now the basic phenomenology is 
rather well-known: noninteracting electrons confined to 
a plane, when subjected to a fixed uniform magnetic field 
and a random scalar potential, exhibit quantum critical 
behavior at a sequence of energy eigenvalues, where the 
correlation length behaves as S,{E) ex \E — E^\~^ with 
i^ ~ I HJll- Early theories of localization in the pres- 
ence of a magnetic field held that no extended states 
exist in two dimensions, a result derived from the cor- 
responding "unitary class" nonlinear sigma model de- 
scription of the long-wavelength physics. Experiments, 
however, unambiguously demonstrated the existence of 
extended states in the quantum Hall system. The work 
of Levine, Libby, and Pruisken [gj showed how a novel 
topological term present in the sigma model could pro- 
duce the necessary critical behavior, although technical 
difficulties rendered the improved sigma model ineffective 
in providing a quantitative description of the transition 
{e.g. critical exponents, scaling functions, etc.) Nonethe- 
less, the language of critical phenomena provided a very 
useful framework within which one could interpret ex- 
periments HH|. Consider, for instance, the behavior of 
the Hall conductivity axy as a function of the variables 
B (magnetic field), L (linear system dimension), and 
(temperature). Within the scaling regime, one can write 

axy{B, L, 6) = ^f{\B - B*\L'/'', \B - B*\ 6-1/^") 

(1) 

where B* is a critical magnetic field, z is the dynamic 
critical exponent, and f{x,y) is a universal scaling func- 
tion. From this expression, one finds the maximum slope 



{d(7xy/dB)^^^ is proportional to 8^^/^"^; the power law 
divergence with decreasing O is in agreement with exper- 
iments, which independently show z — 1 [|6|. Numerical 
calculations JTl confirm the result ly = 2.35 ± 0.02, in 
tantalizing agreement with a compelling but alas nonrig- 
orous argument that i' should be X |q|. 

One of the most significant developments in modeling 
the QHE transition has been the advent of the Chalker- 
Coddington network model of quantum percolation [p|. 
The relevance of classical percolation to the quantum 
Hall problem was emphasized by Trugman ||T^. Clas- 
sical electrons in a strong magnetic field B = Bz obey 
the guiding center drift equations of motion 



dr c 

— = -zx VV"(r) 

dt eB ^ ' 



(2) 



where V(r) is the external potential. A vivid picture 
emerges in which the Fermi 'sea' is analogous to a real 
sea covering a rough surface in energy space whose height 
is described by the function V{r). The corresponding 
quantum eigenstates are localized along equipotentials, 
accruing an Aharonov-Bohm phase 4- ^ A. ■ dv which for 
a complete orbit is proportional to the magnetic flux en- 
circled (single valuedness then leads to Bohr-Sommerfeld 
quantization rules). Electrons at the Fermi level either 
circulate clockwise (as viewed from 'above') around iso- 
lated lakes, when E^ < {V), or counterclockwise around 
isolated islands, when Ep > {V). Such states clearly are 
localized. Only when Ep = E^ = {V) are states at the 
Fermi level extended. As E ^ E^, the equipotentials 
become more and more rarefied, and their circumference 
increases in size as \E — E^\~'^^^, an exact result. How- 
ever, the quantum eigenfunctions are not infinitely nar- 
row. Rather, they have a width on the order of the mag- 
netic length £ = yJhc/eB, and thus quantum tunneling 
will occur in the vicinity of saddle points of V^(r) p,[ll[. 
It is precisely this physics which is captured by the net- 
work model. Thus, rather than directly computing the 
eigenfunctions of 



7i^J-(P + fAf + y(r), 



(3) 



many of which are tightly locahzed on peaks and vaUeys 
and hence irrelevant to the physics of the QHE transition, 
the network model "cuts to the chase" and simulates a 
network of saddle points (see fig. [|), each described by 
an S'-matrix 



O 
O' 



= S 



S = 



t r' 
r t' 



(4) 



relating incoming to outgoing flux amplitudes B . In the 
frame of an incoming electron, scattering is either to the 
left, with probabihty T = \t\^ = |f'p, or to the right, 
with probability R = \r\^ = \r'\'^ = (1-T) - there is 
neither 'forward' nor 'backward' scattering. 

When T < i, scattering is predominantly to the right, 
corresponding to the aforementioned clockwise motion 
around isolated lakes. Hence there is a correspondence 
between the energy eigenvalue E for the Hamiltonian of 
eq. (0) and the transmission probability T; the quantum 
critical point should lie at T = ^. Randomness enters 
the network model principally through phases acquired 
by the flux amplitudes in the course of their propaga- 
tion along each link. These phases reflect the Aharonov- 
Bohm phases accrued due to motion along cquipotential 
segments of irregularly varying length. The phases on 
the links are therefore modeled as random variables uni- 
formly distributed between and 27r. Several numerical 
investigations 0,g,|l2|,Q have convincingly demonstrated 
the applicability of the network model to the quantum 
Hall transition. 



results will be of dubious quantitative value. However, 
we find that a simple and appealing qualitative picture 
emerges from this approach, and we feel its pedagogical 
value alone merits publication. We also wish to draw 
the reader's attention to the recent work of Galstyan and 
Raikh ||lj] , who independently developed a real space RG 
approach to the network model. Their results, of which 
we became aware when this paper was in its final stages 
of preparation, are largely complementary to those pre- 
sented here. 



II. THE DISTRIBUTION Pl{T) 

We commented above how a Chalker-Coddington net- 
work composed of identical scatterers should exhibit a 
quantum critical point when T = R~ ^ . Consider now a 
network in which the individual scattering probabilities T 
(as well as the phases on the links) are random variables, 
chosen from a distribution P{T). We shall be interested 
in how the distribution Pj^{T) behaves with increasing 
cell size L in the limit L —^ oo. To define the transmission 
coefficient T of a finite network, we cut out a \/2L x ^/2L 
section from an infinite network (the lattice spacing is 
a = 1), and retain only the central V{L) = L'^ + {L — 1)^ 
scattering vertices. Fig. g shows the L = 3 case. Note 
that there are M = 2L — 1 scattering nodes across each of 
the main diagonals. Integrating over the random phases 
on the links, one obtains the distribution Pj^iT) for trans- 
mission to the left through the supercell. In the limit 




FIG. 1. The Chalker-Coddington network model of quan- 
tum percolation. 

In this paper, we adopt a real space renormalization 
group (RG) scheme which will allow us to compute ex- 
ponents and scaling functions associated with the quan- 
tum Hall transition. Inasmuch as real space renormal- 
ization is fraught with uncontrolled approximations, our 




FIG. 2. Transmission and reflection through a 3\/2 x 3\/2 
network (L = 3). The 13 vertices are denoted by black dots. 

L — > cx), we expect two stable distributions, given by 
P^[T) = 5[T) and P^{T) = 5[l-T). These correspond 
to bulk localized phases with T — Q and T = 1, respec- 
tively. The quantum critical point will be characterized 
by an unstable distribution Paa{T) — P*{T). The terms 
"stable" and "unstable" refer to renormalization group 



flows ||ij] : we shall develop an approximation scheme by 
which the distribution Pi,]^(T) for a larger system may 

be determined in terms of Pj^{T). This functional rela- 
tion may be represented in terms of a set of parameters 
{Xi{Ly\ which characterize the distribution P]^{T) {e.g. 
the coefficients in a Legendre or Chebyshev polynomial 
expansion of P]^{T) in the variable x ~ 2T — 1): 



X,{bL) ^ F,{{X,{L)};b) 



(5) 



When X* — Fi{{X*};b) the distribution is at a fixed 
point. The eigenvalues {Afe} of the matrix Rij — 
{dFi/dXj)\^^ determine the relevance of the correspond- 
ing eigenvectors (scaling variables) and hence the stabil- 
ity of the fixed point - positive eigenvalues correspond to 
relevant variables and an unstable fixed point, negative 
eigenvalues to irrelevant variables. The positive eigenval- 
ues define a set of critical exponents, y^T = In 6/ In A^; 
the localization length exponent v = Inb/ In Amax is the 
largest of these. One can also define the beta function 



mxj}) - 



dX, _ dF, 
d\nL db 



6=1 



(6) 



the beta function vanishes at all fixed points. 

From these considerations it becomes clear that one 
possible way to compute the critical exponent v would be 
to determine the RG flow of P^ (T) . An exact procedure 
would require the calculation of this distribution for finite 
networks of arbitrary size. Computational limitations 
would render such attempts futile beyond even modest 
values of L. Instead, we shall resort to two uncontrolled 
and closely-related approximations: (i) Migdal-Kadanoff 
(MK) bond shifting, and (ii) replacing the square lat- 
tice with an unphysical, hierarchically constructed lat- 
tice. These approaches will allow us to perform the renor- 
malization recursively, using only simple numerical com- 
putations. 



III. SERIES AND PARALLEL PROPAGATION 

The group U(2) of unitary 2x2 matrices may be pa- 
rameterized by four angles (a, f3,j,0): 







= -»/3 



COS 9 sin t 
— sin cos ( 



g«7 





-27 



The phase angles a, /3, and 7 can be absorbed by the 
random link phases, and hence without loss of generality 
one may restrict attention to 5'-matrices of the form 



S = 



cos 9 
— sin 9 



sm t 
cos( 



(7) 



This S'-matrix may be recast in the form of a "left-to- 
right" transfer matrix Af , which relates flux amplitudes 
O and T to I and O': 



I —^ 




0'^ 



-^0 



(b) 



I 



0' I' 

(c) 



FIG. 3. (a) The S-matrix relates incoming to outgoing 
flux amplitudes (see eq. (ft)). This information may alter- 
nately be expressed as a transfer matrix M relating left and 
right flux amplitudes (b), or A'^ relating top and bottom flux 
amplitudes (c). 



= M 



A/ = 



I 
O' 

sec 9 tan 9 
tan 9 sec 9 



The transmission coefficient T relating the incident flux 
|Xp to the outgoing flux \0\^ is T = cos^ 6*, and the 
reflection coefficient is i? = 1 — T = sin^ 9. 

Consider now the "series" transmission through two 
consecutive scatterers. Multiplying transfer matrices, 
and taking into account the random phases accrued in 
between scatterers, one obtains a composite transfer ma- 
trix 



Af = 



' sec 6*2 
tan 6*2 



tan 6*2 
sec 02 











' sec 9\ 
tan 6*1 



tan I 
sec( 



and from scc^ 9' = |A/i,ip = 1/T', a transmission coeffi- 
cient 



1 



1 1 



1 + 2 v/i?ii?2 cos (j) + RiR2 



(8) 



where (/> = a — /? is random. Precisely this calculation 
was done by Anderson et al. pq ] in their study of one- 
dimensional localization (see also refs. [^ 18|. Averaging 
InT over the angle (/>, one obtains 



(InT') =b{\nT) 



(9) 
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FIG. 4. The horizontal 'bricklayer's representation' of the 
network model. Quantum mechanical particles propagate in 
the directions of the arrows; scattering processes are repre- 
sented by the dotted lines. The corresponding vertical repre- 
sentation is rotated by 90°. 



Thus, InT is driven to in- 
- this is the essence of one- 



for b scatterers in series, 
creasingly negative values 
dimensional localization. 

Equivalently, though, one may view the propagation as 
"top-to-bottom" , and define a transfer matrix N relating 
C and T to X and O: 



O' 
I' 



N 



N : 



o 

- CSC 6 ctn 9 
-ctn6' CSC 9 



In this case, the roles of T and R are reversed: T 
R — cos^ 9, and two scatterers in parallel give 



1 
R' 



1 
R'i 



1 



with 



l + 2^/TiT2COS(j) + TiT2 



(Ini?') =5(lni?) 



(10) 



for b parallel scatters. In this case the reflection am- 
plitude is driven to zero. Now in the network model, 
both series as well as parallel propagation occur, and 
in our viewpoint it is a competition between these two 
one-dimensional localization mechanisms which leads to 
a critical point corresponding to the quantum Hall tran- 
sition. It is, however, impossible to neatly separate the 
two modes of propagation, and we shall have to resort 
to approximation schemes, such as the Migdal-Kadanoff 
approach, or to a modification of the original model, such 
as hierarchical network constructions, in order to make 
progress. 



FIG. 5. The first step in the 6 = 2 Migdal-Kadanoff 
renormalization of the network model. Compare with fig. H. 
Dashed lines indicate bonds which have been shifted to the 
left. Dotted lines indicate unshifted bonds. 



IV. MIGDAL-KADANOFF APPROACH 

In order to apply the Migdal-Kadanoff bond shifting 
scheme to the network model, we begin by graphically 
representing the elementary scattering process as a kind 
of interaction vertex, as depicted in fig. |[ Thus, the en- 
tire network may be represented by an infinite sequence 
of parallel lines on which currents flow in alternating di- 
rections, and backscattering processes which divert flux 
from one line to one of its neighbors (see fig. 0). We 
call this the 'bricklayer's representation' of the network 
model, for obvious reasons, and further describe the prop- 
agation as being either 'horizontal' or 'vertical'. To effect 
a rescaling by a factor &, we start with the (horizontal) 
bricklayer's representation, shown in fig. H. In every other 
row, we then shift 6—1 out of every b interaction lines so 
as to form a new network in which each renormalized ver- 
tex represents b bare vertices in series. This is depicted 
for 6 = 2 in fig. H. We then rescale horizontal distances 
by b. Were we to proceed in this manner, we would ob- 
tain the one-dimensional localization of rcf. pq. Instead, 
we view the network again in the bricklayer's representa- 
tion, this time with vertical propagation, and perform a 
second bond-shifting operation. This replaces each scat- 
tering vertex with an effective composite vertex, shown 
in fig. 0. The new network resembles that of fig. n 

To see how this leads to critical behavior, consider the 
behavior of the typica? transmission coefficient, exp(lnr), 
under such a transformation. From eq. (y), we obtain 



T' = 1- (l-T 



b\b 



/(T;6) 



(11) 



where we implicitly are working with typical values. For 





FIG. 6. The composite vertex in the b — 2 

Migdal-Kadanoff renormalization, which corresponds to se- 
ries followed by parallel propagation. 



FIG. 7. After shifting in both the horizontal and vertical 
directions, each vertex in the network represents a composite 
of b^ vertices. The case & = 2 is shown here. 



the case 6 = 2, we have 

/(T; b = 2) = 2T^ -T'^ , 

which has two stable fixed points at T* = and T* = I, 
and one unstable fixed point at T* — i(\/5 — 1) ~ 0.618. 
Linearizing about the unstable fixed point, we obtain an 
eigenvalue {df{T; b = 2)/dT)j,, = 6 — 2\/5, correspond- 
ing to a localization length exponent oi v — Inb / \n\ k 
1.635. Note that T* 7^ i, a consequence of the order 
in which the bond shifting was performed: the com- 
posite vertex corresponds to series followed by paral- 
lel propagation (note the composite vertex of fig. || are 
not symmetric under 90° rotations). We could equally 
well have chosen parallel followed by series propagation, 
which would switch the roles of T and R. We note that 
the RG equation (O) and its counterpart with T replaced 
by 1 — r coincide with two RG equations obtained in a 
Mig dal-Kadanoff approach to classical bond percolation 
|l9| . There, the bond occupation probability p plays the 
role of T, and the MK bond shifting which leads to the 
composition of quantum resistors in the network model 
corresponds to multiplication of bond occupation proba- 
bilities. 

If we set b — \ + C,, where ^ — *■ 0, we obtain equations 
for an 'infinitesimal' Migdal-Kadanoff renormalization: 

T' = T + C[ThiT - (l-T)ln(l-T)] + C'(C^) . (12) 

The infinitesimal MK transformation has a fixed point at 
T* = i, and an eigenvalue A = 2(1— In 2)^, corresponding 
to a critical exponent of i' = 1/2(1 — In 2) « 1.629. The 
MK beta function is then 



dC, 



/(T;l + C) 



C=o 



= rinT- (l-r)ln(l-T) 



(13) 



this is shown in fig. g|. 

In statistical mechanics applications, such as the Ising 
model, one can define a new Hamiltonian TC in which 
bonds have been shifted according to the MK prescrip- 
tion. When computing the partition function, one can in 
principle perform thermodynamic perturbation theory in 
(H' — H) in order to systematically improve upon the MK 
procedure. We know of no such systematic improvement 
for our scheme, nor do we have any sort of reliable esti- 
mate for the errors involved in the calculation of v. 



V. HIERARCHICAL LATTICES 

A related approach to the problem involves the con- 
struction of hierarchical lattices. Consider, for example, 
the scattering unit of fig. ^ which contains V = 5 vertices 
in an area M^; in units of the distance between vertices, 
we have M = 3. Now replace each of the vertices with a 
replica of the original cell, forming the structure shown 
in fig. |o[ Repeating the process n times generates a hier- 
archical structure with V"' vertices contained in a square 
of side length i\f " . The Hausdorff dimension of this hi- 
erarchical lattice is I?jj = \nV/\nAI « 1.465. From 
ViM) = i(M2 + 1), one has limM^oo D^iM) = 2. The 
renormalization flow of the transmission coefficient and 
the multifractal spectrum /(a) was calculated in ref. pQ] 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 



T 

FIG. 8. The /3-function resulting from the infinitesimal 
Migdal-Kadanoff renormalization. T = and T = 1 are at- 
tractive fixed points; the critical point lies at T = ^. Arrow- 
heads point in the direction of RG flow. 



for the cases M = i, M — 5, and M — 9; the resuhs for v 
and f{a) rapidly converge to the network model results 
even for such modest sizes. The RG equations for the 
case M — i are identical to those obtained in ref. p4[ , 
where a different interpretation is adopted. 

If one replaces the central scatterer in fig. |^ with one 
for which T = 1 , one recovers the four site scattering unit 
of fig. o. One can use this as the fundamental unit of a hi- 
erarchical construction, and the results differ from those 
of the previous section only in that the linear dimension 
is taken to be M = 26 — 1 rather than b. The Hausdorff 
dimension of the hierarchical lattice is D^ — Info'^/lnM 
{— In 4/ In 3 for 6 = 2), whereas the lattice of fig. is 
fully two-dimensional. The exponent u now is changed: 
J/ = InM/lnA = ln3/ln(6 - 2^5) « 2.592. Generaliza- 
tions to other values of h are straightforward: iterating 
eq. (O) n times, we obtain the flow equation 



r(M"4)=/["l(r(4);6=iM 



(14) 



where /'"' is the iterated function satisfying 



/["+il(r;6) = /(/N(r;6) 



and /[^l = /. Since /(T; 6) is monotonically increasing 
from /(0;6) = to /(I; 6) = 1, there are always three 
fixed points on the interval T G [0, 1]: T = and T = 1 
(both stable), and a nontrivial unstable fixed point with 
< T* < 1. The critical exponent v is 



InM 



ln/'(T*;iM+i) 



(15) 
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FIG. 9. An M = Z,V 
hierarchical lattice is based. 
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5 elementary cell on which a 



where /'(T; 6) = df{T; b)/dT. 

Whereas the Hausdorff dimension D^ 
both increase monotonically, with 



and fixed point 



D^ib 

r*(6 = 



D^{b 
T*{b^ 



oo) = 1 



the critical exponent v exhibits a minimum for 6 « 35, 
where i' « 2.11. We find that i' diverges weakly (as 
In 6/ In In b) in the b ^ oo limit. These results are plotted 
in fig. 0. 



VI. 



RENORMALIZATION GROUP FOR 
DISTRIBUTIONS 



In section IV , we derived RG flows for the typical trans- 
mission coefficient, defined by Ttyp = exp(lnr). Now we 
will concern ourselves with the RG flow for the entire dis- 
tribution P{T). We will choose the simplest nontrivial 
unit cell, namely that of fig. 0, in which two composite 
scatterers, each of which is two bare scatterers in series, 
are placed in parallel. 

Starting with a distribution function P{T), we derive 
the intermediate distribution Pj^^^ after taking two scat- 
terers in series: 

^int(^) = fdT.PiT,) fdT.PiT^) X 
Jo Jo 

X / ^J(T-/(Ti,T2,cos0)) (16) 

Jo "^ 



-^ — II — f- 
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FIG. 10. The hierarchical lattice at the second level of 
construction, beginning with the M = 3, V — 5 cell of fig. bl 
There are V^ = 25 vertices in this M^ x M'^ =9x9 supercell. 



FIG. 11. Results of the hierarchical lattice renormahza- 
tion of eq. (|l4) for the fixed point transmission coefficient, T* 
(solid curve), and the localization length exponent, v (dashed 
curve), versus In 6. The linear dimension is M = 2& — 1. 



with 



I{T,,T2,y) 



T1T2 



1 + (1 - Tl)(l - T2) + 2^1y/il-Tl){l-T2) 



We now combine two of these composite units in par- 
allel, obtaining the renormalized transmission coefficient 
distribution 

P{T) = [ dR, Pi„,(l - R,) fdR2 Pi„,(l - R2) X 



X / -1^(1 -T-/(i?i,i?2, cos 0)) . (17) 

Jo 71" 



This is equivalent to the following: 
PintiT) - ^[dT^J^dT^ (Ti) P{T2) 



/o Jo 

T1T2 



e(i-|A|) 



A = 



v/(i-ri)(i-T2) 
TiT2-r(i + (i-ri)(i-r2)) 

2Ty/{l-Tl){l-T2) 



followed by 



^^ 9(i-|i3|) (i-ri)(i-r2) 



B = 



(1 - ri)(i - T2) - (1 - r)(i + Ti r2) 
2(i-r)V7YT^ 



Note that B{Ti,T2;T) = A(l - Ti, 1 - T2; 1 - T). 

We have numerically iterated the distribution accord- 
ing to these equations. Since we expect a single rele- 
vant scaling variable controlling the RG flow of P{T), 
we parameterize the initial distribution by its average 
(r)i„i, and a width (AT).„., = ^{T^-{Ty). By vary- 
ing (r)jjjjj alone, we drive the system through its critical 
point. This is shown in fig. [ij, where the iterated aver- 
age transmission coefficient is plotted for several initial 
values of (T). (We used an initial distribution which 
was a sum of linear plus exponentially increasing terms, 
though the precise details are presumably unimportant. 
Keeping the initial width fixed and varying (T')jj^j^ be- 
tween 0.605 and 0.606, we found the critical point (T)* 
at (T)jj^jj Ri 0.60536.) The fixed point distribution itself 
is plotted in fig. Ol 

Within the hierarchical lattice approach, each iteration 
rescales the linear system size by a factor M = 3. An- 
alyzing the iterated average transmission coefficient (T) 
as a function of the number A'' of iterations allows for the 
computation of v from 



(T)^-(r)*(x Af^/" . 



(18) 



In fig.|l4|we plot In |((r)jv - (T)*)/(r)*| versus NlnM; 
we find zy « 6.6. Note that within the MK scheme, the 
rescaling factor is & = 2, consequently v « 4.2. 




ITERATION 

FIG. 12. The RG flow of the distribution P{T). The 
average transmission coefiicient (T) is plotted versus iteration 
for various initial conditions. The initial value of AT was 
held flxed as {T}-^^^^ was varied between 0.605 and 0.606. The 
critical point was found to lie at {r)j^jj « 0.60536. 



The iteration of the full distribution function P{T) 
according to eqs. ( p^ , p7[ ) yields an exact RG flow for 
Pi^iT) for the four site scattering unit of fig. ^. How- 
ever, the iteration could be performed only numerically. 
We conclude this section with a brief account of an an- 
alytical method to obtain P{T). In the context of one- 
dimensional localization eq. (pf) has been used to derive 
a Fokker-Planck equation (c/. [gl[) for the distribution 
function P^{T). One identifies Ti = T{L) with the trans- 
mission coefficient of a one-dimensional resistor of length 
L placed in series with a resistor of infinitesimal length 
and transmission coefficient T2 = T{SL). In a regime of 
'local weak scattering', one can identify the elastic mean 
free path l^ by le = SL/{R{SL)) and derive the follow- 
ing Fokker-Planck equation by standard methods (c/. p^ ] 
and references therein): 



° OL drV ' dT 



T^Pl{T) 



(19) 



This is an equation describing one-dimensional localiza- 
tion. It turns out that Ic sets the scale for the localization 
length, i.e. (-InT)^ = L/l^. 

We now try to extend this one-dimensional scheme 
to the network model. Our two-dimensional iteration 
scheme is graphically depicted in fig. HS. We begin with 
an L X L array and attach to it two infinitesimal rectan- 
gles of size LxSL and SLx (L + SL), respectively. Trans- 
port then takes place according to the usual series and 



[^ 




FIG. 13. The fixed distribution P*{T) obtained by it- 
eration of the RG flow equations. Shown are results of a 
numerical integration (dotted curve) and a smooth fit to the 
numerical data (solid curve). 



parallel composition laws. The main obstacle in deriving 
a Fokker-Planck equation for this model is due to the 
fact that we do not know a priori how the transmission 
coefficients of the infinitesimal blocks behave. To make 
progress, we proceed in analogy to the one-dimensional 
case and assume 



{R{L,SL)) = {-\nT)^-{SL/L) 
(T{SL, L + 5L)) = {- ln(l - T))^ • (SL/L) 



(20) 



With this assumption the resulting Fokker-Planck equa- 
tion reads 



d 



d 



^^.(r)^--:.(r;x). 



where the 'current' is given by 
l-T d 



^(r;L) = (lnr)^. 



2i dT 



{T'P^iT)}] 



(21) 



(22) 



T d 



+ {ln{l-T)),- — {il-TrP,{T)}. 



Without going into details we report on the results that 
one can obtain from the Fokker-Planck equation (pTl) . We 
find that there are three fixed point distributions: S{T), 
(5(1 — T), and a uniform critical distribution P*{T) = 1. 
To extract the critical exponent v of the localization 
length we took Xj^ = (ln(l — T))^ — (InT)^ as a scaling 
variable which has a fixed point value X* = 0. Lineariz- 
ing the flow equation around the fixed point distribution 
gives rise to a beta function for Xj^ which reads 



z 





Nln3 



L + 25L 



FIG. 14. To derive v from eq. hq, we plot the logarithm 
of the left hand side versus NlnM. The slope is then v~^ . 
Here we show results for M = 3 discussed in sec. W% we find 
v « 6.6, in poor agreement with the accepted value u ^ ^. 






(23) 



Consequently, the critical exponent is i^ = 1. Although 
the Fokker-Planck approach to the network model is 
able to describe the correct qualitative physics of the 
localization-delocalization transition, it suffers from an 
ambiguity in modeling the average transmission coef- 
ficients of the infinitesimal blocks. Adopting a dif- 
ferent model can change the critical exponent v. In 
addition, the composition laws on which the Fokker- 
Planck approach is based do not lead to the full Chalker- 
Coddington network. 



VII. CONCLUSION 

In this paper we developed a simple approach to 
the quantum Hall transition. We considered a number 
of hierarchical network models related to the Chalker- 
Coddington model of quantum percolation. The basic 
ingredients of our networks are series and parallel com- 
positions of quantum resistors; their hierarchical nature 
allows us to analyze their critical properties using simple 
real space renormalization techniques. In particular we 
studied network models that can be interpreted as result- 
ing either from Migdal-Kadanoff bond-shifting scheme or 
from the construction of hierarchical lattices of Hausdorff 
dimension D^ < 2. The different models are labeled by 
an integer b where b^ is the number of resistors that form 
an elementary cell. 

In sec. pll we calculated the flow of the distribution 
function P^T) for the transmission coefficient T ('two- 
point-conductance') in the case 6 = 2. This was done by 
means of numerically iterating the renormalization group 



FIG. 15. Construction of a network by adding infinites- 
imal networks to a square network of linear size L. The in- 
finitesimal networks are characterized by a small extension 5L 
in one direction. The network is composed according to series 
and parallel composition laws. 



equations. We found the flow to have three fixed points, 
two of which correspond to localization (T = 0, T = 1). 
The third fixed distribution P*{T) corresponds to the 
quantum critical point, where delocalization occurs. The 
critical point distribution P*{T) turned out to be very 
broad and the flow of P^ (T) in the vicinity of P* (T) is 
governed by one-parameter scaling. The critical expo- 
nent V was also determined and found in poor agreement 
with accepted values. 

As a further simplification we derived renormalization 
schemes for the typical transmission coefficient, Ttyp = 
exp(lnT). Here, we simply averaged over the random 
link phases, neglecting variations in the individual trans- 
mission coefficients. We obtained closed RG equations for 
Ttyp (sees. IV, ^. From these equations we derived fixed 
point values for T^ (b) and critical localization length 



exponents i^f, for arbitrary & > 1. For small b the no- 
tion of typical transmission coefficient as a substitute for 
a whole distribution (which is broad) is dubious and it 
is no surprise that i'j^2 differs from the value obtained 
by iterating the whole distribution function. However, 
the notion of typical transmission coefficient is more reli- 
able for large b where the composition of b resistors in se- 
ries favors the formation of log-normal type distributions. 
On the other hand, for large b the resulting hierarchical 
network has little in common with the original Chalker- 
Coddington network. Nevertheless, the slow variation of 
Vj^ with b shows that the results for intermediate b are 
not too far from that of the original network model. 

All models that we considered exhibit a localization- 
delocalization transition which results from a competi- 
tion between two one-dimensional quantum mechanical 
localization mechanisms. The models result from uncon- 



trolled approximations to the network model, but have 
at least some pedagogical value. In contrast to the well- 
known classical percolation model for quantum Hall tran- 
sitions, our models are essentially quantum mechanical. 
Furthermore, the concept of real space renormalization 
can be extended to improve quantitative results (see refs. 

While this manuscript was in the final stage of prepa- 
ration, we discovered the work of ref. |1J| in which a 
similar real space renormalization group approach to the 
Chalker-Coddington model is developed. In that work 
an elementary M — S, V = b cell (see fig. ||) replaces 
each scatterer in the renormalization step. (In addition, 
the rescaling factor in ||lj| is 6 = 2 rather than M — 3.) 
In our work the central scatterer of this cell is replaced 
with one for which T = 1 (or _R = 1), resulting in the 
simple series-parallel composition laws discussed above. 
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